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Dynamic  Response  and  Control  of  Multipass  Heat  Exchangers 


MASAMI  MASUBUCHI 


Introduction 


Pure  parallel-and  counterflow  heat  exchangers  have  been  analyzed 
(1 )  (2^ 

exactly  and  their  controllability  have  been  extensively 

discussed  in  the  literatures  For  multi-pass  heat  exchangeraulv^.  . 

nany  kinds  of  researches  on  the  mean  temperature  difference 

n  1  ^ 

and  on  steady  temperature  distributions  '  '  have  been  published » 

but  there  have  been  few  on  dynamic  analysis* 

However,  from  the  control  engineer's  point  of  view,  process  dynamics 
is  one  of  the  most  Important  subjects. 

The  dynamic  response  of  the  system  is  dependent  not  only  on  the 
mean  temperature  difference,  but  also  on  other  system  parameters  such 
aJ^  fluid  velOT^fties,  heat  capacities  of  the  fluids  and  tube-walls, 
etc.*  The  transient  response  seems  quite  difficult  to  calculate,  but 
the  frequency  response  analysis  is  simpler  and  is  used  throu^out 
this  paper. 

The  author  extended  the  P, Prof os's  basic  partial  differential 
(ip) 

equation  '  '  and  analyzed  exactly  the  dynamic  characteristics  of 

1-2, 3, 4,  -  ,2n,2ttH  pass  heat  exchangers. 

The  3rd  order  characteristic  equations  goveming  the  heat  exchange 
processes  in  each  multi-pass  heat  exchanger,  and  the  transfer  fimction 
formulas  which  involve  the  roots  of  the  said  equation  are  obtained 
and  presented. 

Numerical  examples  of  frequency  response  of  each  case  with  no  solid 
capacity  have  also  been  shown. 


Nomenclature 

The  following  nomenclature  is  used  in  this  paper: 

2 

m  *  heat  exchange  surface  area 


kA/q,  a^»  kA/q^^ 


^Numbers  in  parentheses  refer  to  similarly  numbered  references 
in  bibliography  at  the  end  of  the  paper. 


2 


3 


/ 

a 


dl. 
b  dl. 

C  lccal/n‘’C 
H  m 

h  kcal/m^mn‘*C 


=  hA/q,  a^=  h^^A/q^,  3^=  ligA/q 

=  hA/(Gv),  (CjjV^ ) ,  b2=h]^A/ (Cj^v^ ) ,  bg*hgA/(CgV^) 

=  series  and  side  capacity 
»  total  length  of  shell 

»  film  coefficient  of  heat  transfer,  h^  tube-side, 

h  shell-side,  h_  side  capacity-side 
s 


2 

keal/m  mn^C  *  overall  coefficient  of  heat  transfer 
mn  *  dead  time,  H/v»L,  H/Vj^= 

dl.  »  Y/H 

«  variable  in  axixiliary  equation 


kcal/mn"C 

dl. 

dl. 


or  m/an 
w  kcal/n‘‘C 


o^,  a2,  (Xj 
^ »  ^2  *  ^5 

0  -c 

e  ®c 


'  qi* 

=  Vj^/v  =  L/L^ 

*  jw  ,  j  * 

*  time 

«  peripheral  length  of  heat  exchange  siurface 
e  variable  (  see  Appendix  ) 

«  variable  (  see  Appendix  ) 

*  fl\iid  velocity,  v^^  tube-side,  v  shell-side 
»  heat  capacity  of  fluid  per  tinit  length 

*  ruxming  length  of  shell  measured  from  shell  input 
»  parameters  defined  by  Equation  (3) 

«  parameters  defined  by  Eqxiation  (3) 
s  temperature  of  shell  fluid  in  lumped  system 

*  temperature  of  fluid,  ©  shell- aide,  0^,02,©^, - 

tube- side 


*0  *=  tea^erattire  of  solid, ^  series  capacity,^  side  capacity 

dl.  *  circvilar  frequency 
and  dl .  dimensionless 

Subscripts:  l,2,3f  —  ®  tube-side,  no  subscript  =  shell-side 

1  =  input,  0  *  output 
h  =  series  capacity,  s  *  side  capacity 

Basic  Assumptions 

1.  Flxiid  velocities,  heat  transfer  coefficients  are  all  constant  and 

do  not  change  with  the  temperature  of  fluid  or  heat  exchange  surfaces. 

2.  Each  cross  section  is  constant. 

3.  Complete  mixing  in  crosswise  directions  of  each  flow  exists. 
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□  a 


/ 


4.  The  heat  conductivi  of  tube^-and  shell-walls  are  infinite  in 
the  directions  at  ri  angles  4o  the  flow  and  zero  in  flow 
directions. 

5.  There  is  no  heat  loss  to  outer  circumferences. 

The  transfer  function  is  taken  to  be  the  dynamic  ratio  of  the 
outlet  temperature  of  one  fluid  to  the  inlet  temperature  of  another 
fluid.  Considering  sinusoidal  solution,  sinusoidal  temperature 
variation  is  superposed  on  the  steady  state  temperature,  so  the 
above  definition  is  enough  to  consi'^r. 


Fundamental  Equations 


The  following  relation  exists  between  a  fluid  at  temperature  6 
and  a  solid  surface  at  temperature  ^  through  which  heat  exchange 
takes  place. 


M 


Using  nondimensional  £  and  af  ,  gives 

If  the  tube-wall  is  infinitely  thin,  can  be  replaced  by  the 
temperature  of  another  fluid,  and  replaced  by  a.  This  fundamental 
equation  applies  to  - in  the  following, 

1-2  Pass  Heat  Exchanger 


Case  P-C: 

The  "P-C”  flow  pattern  is  shown  in  Fig,l,  In  this  ease,  the  tube 
fluid  flows  first  from  left  to  right,  turns  aroimd  at  the  right  end, 
and  finally  flows  back  to  the  left,  resulting  in  a  parallel-flow 
followed  by  a  counter  flow.  The  shell  fluid  is  always  assTined  to 
flow  from  the  left  end  to  the  right  end.  Both  inside  and  outside  tube- 
walls  are  assumed  to  be  infinitely  thin  so  that  there  will  be  neither 
series  nor  side  capacity  to  be  considered. 

Asaruaing  equal  dimensions  and  conditions  for  each  tube-pass  axid 
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P  c 


«11®‘  ee,  0,  SHELL  SIDE  FLUID 


'  BAFFLE  PLATEV 

Plg.l  1-2  pass  heat  exchanger  snow¬ 
ing  symbolism. 


0,_o  e. 

(Q)  (b) 

Fig. 2  Case  P-C  (a)  shell-side 
Input;  (b)  tube-side  Input. 


equal  TJ  for  shell-and  tube-fluid,  the  fundamental  equations  are, 

'  ^  +  If 

■  II 0 

,r||+||=a(0,+e)fa(era) 

The  solution  of  Equations  (1)  for  a  sinusoidal  input  signal  (0^) 

•  e,=  (oi|C,^+o(4e^ +0(346^)6^^ 

where,  C^,  Cgt  are  the  integration  constants,  and  a,  ^  are: 


ol  r- - -  qI  =  -r — ^ 

-H  +  stJi,  '  -t>x  +  s 


•vs  -vai 


fe+s  ti, 


Z-i  a  -  w _  Q _ Zb _ 

f»i-s+a,,'  pj+s  +  a,  >  ft' f>3+-s+a, 


P2»  P3  three  roots  of 

t^  +  K(2a+Ts)-Ks+a.T-«U  +  a,){t(s+ai)+2a}=  0  (4) 

The  roots  are  all  complex  (  see  Appendix  ) . 

Shell-side  input:  The  C^,  C^,  can  be  determined  from  the  following 
boundary  conditions  (  Pig. 2(a)), 

Ji  =  0  :  ^  =  Pi(i»  +Pz<^2  -  0 

1=  1 ;  6,  =  e,,  or^,C,e*'  +  |^cl,e^+|^4e^ 

=  of,c,e**  +  ofj^efe 

The  substitution  of  the  results  thus  obtained  into  Equations  (2) 
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9ii“0  0i~O 


a.  %  0  9'  0' 


9k,  0  0  0 


IQ)  (Q)  Pig. 4  With  side  capacity. 

Pig. 3  Case  C-P  (a)  shell-side  in¬ 
put;  (b)  tube-side  Input. 


%  % 

Pig. 5  With  series  and  side 
capacity . 


and  solving  for  ®2o/®i  gives  the  transfer  function, 


&)e^+(o(r 

•fetXfej 

+ (di' 

If  input  and  output  are  exchanged  with  each  other  (Fig.2(b)),  we 
obtain  the  following  transfer  function  through  a  similar  procedure : 


(’06-K<j8>--A)e  ^ 

(I 

Static  Characteristics: 

This  is  the  case  when  s=0  in  Equations  (5)»(6),  For  s=0,  from 
Equation  (4),  three  roots  are 

which  give  ^  a 

Introducing  these  values  into  Equations  (5)»(6),  and  rearranging, 
we  obtain 


_  o,  s  ■ 


^  ^  where, 

(5)»(6)  are  the  static  characteristics  known  as  temperature  efficiency 

(14) 

of  the  heat  exchanger  ' 


Case  C~P: 

Counterflow  at  input  aide,  parallel-flow  at  output  side  (  See  Fig, 3  ). 
The  fundamental  equations  are  as  follows, 

•  ^  =a,,(&-Q) 

.41  +||=al6,-e)+aM 

By  solving  these  eqviations,  the  3rd  order  characteristic  equation 
identical  to  Equation  (4)  can  be  obtained.  The  temperatures  are: 

'  d,  =  (a^;e^+ctC;^e'‘‘^+oi^e^*-) 

where,  C^,  €2$  C^  are  the  constants  of  integration,  a,  ^  can  be 
expressed  by  the  same  equation  as  (3)*  Through  the  same  procedure  as 
before,  we  get 

(9) 

KCp.-  «»-  <fe-cfe)>«+ 

The  static  characteristics  are, 
fro.  Equation  (9)  C^cp-sttO) 

from  Equation  (10)  _  ^^/n\ _ _ 

^cF-iscoj-  a,+a.tvff'c«thvg^ 

With  Side  Capacity 

Assume  that  the  side  capacity  (  thermal  capacity  of  shell  wall  ) 
has  infinite  thermal  resistance  in  the  axial  and  no  thermal  resistance 
in  the  radial  directions. 

In  the  case  P-C  (  same  as  C-P  ),  for  example,  the  fundamental 
equations  are  (  See  Pig. 4  ) 
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(11) 


-  Hh-II  =a.(e-0,) 

II  =W(e-!6) 

In  this  case,  the  characteristic  equation  (4)  changes  somewhat 
and  takes  the  following  form 

f  (2  a.  tTs  + -^ )  -  Ks 

-s(s+a,){v'(s+a,)t2a+^(s-(-a,)}  =o  (i2) 


Viith  Series  and  Side  Capacity 

Assxime  that  the  series  {  solid  capacity  of  tube-wall  )  and  side 
capacity  have  infinite  thermal  resistance  in  the  axial  and  no  resistance 
in  the  radial  directions.  In  the  case  of  P-C,  for  example,  the 
fvindamental  equations  are  (  See  Pig. 5  )» 


1^+'^  =^U%-e,) 


(13) 


From  the  above  equations,  the  3rd  order  characteristic  equation 
with  complex  coefficients  can  be  obtained.  The  result  is: 

)?+wirts^.^  +2Aktfe»)i_b  r^^4iTstt>,)  -f 

^  VK+b*)^^S+i,->;b. 

and  Equations  (3)  changes  to 


(14) 
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(t»  (C) 

Pig. 6  Case  P-C,  1-3  pass  heat  ex¬ 
changer. 


e  0,  0/  0> 


Om 


0. 


(b)  (C) 

Pig, 7  Case  C-P,  1-3  pass  heat  ex¬ 
changer  . 


S  +  h  +  tr 


a; 


r.  h 


— a —  2^ 


hs+ 


s±bi 


■af, 


(15) 


€+!)(  +t>t' 

The  transfer  fiinctions  may  be  obtained  through  a  similar  procedure. 


1-3  Pass  Heat  iixchanger 


Case  P-C-P: 

This  is  a  parallel-counter-parallel  flow  heat  exchange  process. 

The  fundamental  equations  are  as  follows,  when  no  solid  capacity  exists, 

'  ^  +  =a.{e-9.} 

,  =i.C0-0,) 

.'"H  +  il  -a(e-e)+!>.{e,-d)+a(9j-d) 

The  temperature  variations  of  each  fluid  when  input  fluid  temperature 
change  is  sinusoidal  (  )  are : 

where,  D^,  D2»  ^3»  are  the  constants  of  integration,  and  P2,P2»P3 
are  the  roots  of 
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l^  +  ^(3a.+rs)-|>{(Ji,+sf +WtJ-S(s+AO(T(^,+  s)  +3a}=0 

and  ‘  ^  ^  1 

The  a,  (3  are  given  by  Equation  (3). 

If  shell  input  is  considered,  from  the  similar  boundary  conditions 
as  1-2  pass  type,  as  shovn  in  Fig.6(b).  the  constants  of  integration 
are  determined  and  the  transfer  function  &pcp-st^®^  "  ®3o/®i 
obtained:  . 

^pcp-st  CS)  =[Waft  ) 

where,  t 

*♦-  e^(0(,-  Bj)[2  -  (0(r<Vr)}  +  e*^[Q(,  (f2“^3)t 

(21) 


Changing  input  and  output  relations  as  in  Fig,6(c),  we  get 

^PC-tsCs-) 


(22) 


Case  C-P-C; 

This  is  the  case  of  cotinter-paralle  1-counter  flow  heat  exchange 
process.  See  Fig, 7(a).  The  fundamental  equations  are: 

'  ^-t=a.(e-e.) 

^  11  +  11=  a, (e-&) 

II -If  =a,(0-Q.) 

+^  +  ||  ca(6!-©)+4.(ez-e)ta(%-e) 


From  which  we  get 

r  0  =  (<v,p,e*''®’+o(iPxe^+oijr^e^+ 

I  k_  A  i  A 


9  =  Aire®V* 

t_  />  1-  A  L  n 


I  0  ~  ( i>.e^^+  Pxe^+P,e^*')e^^ 


(23) 


(24) 
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where,  D2»  are  the  integration  constants  and  p^,  p^,  p^ 

are  the  three  roots  of 


f (3a.+'^s)-i.|(4,+sf-aa,^-s(s+^,){r(a,+s)  +3a}  =  0 

and  ■  a-  +  s 


(25) 

(26) 


5  -  ^  ® 

01,  p  can  be  expressed  by  Sqtiation  (3)  as  before. 

Shell  side  input:  Fro*  the  boundary  conditions  shown  in  Pig.7(b), 
we  gtt  the  following  transfer  fxmction  throvigh  a  similar  procedure 
as  before. 

+  (2<r,-)5,Xa<.^  e*' 

where,  ■ 


(27) 


Tube  side  input:  From  the  boundary  conditions  in  Fig, 7(c),  we  get 
the  following  result: 

(ot-A^e^-e^'^'^lA;  (29) 

1>4  Pass  Heat  Exchanger 

Case  P-C-P-C:  Assuming  as  before,  the  fundamental  equations  are: 

=Md-tk) 

Ic  +  m  ^  1?)  ^A(£v-e) 
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0  0  0,6=0.  0 


Fig. 8  Case  P-C-P-C,  i-4  pass  heat  Fig -9  Case  C-P-C-P,  1-4  pass  heat 

exchang  er .  exchang  er . 


Froai  which,  we  get 

-  e,  =(^E  ■“'+&E.e^‘'tAE3e^“'+Afe^) 

61,  =(o|,E,e^*'+«»Exe^‘^+o(jEie^- 

,  a,  =(eiE,e^^E,€!^+^Eje^- ‘ 
fiV  =(o(,E,«?'’+ofiEie.^'^+«5£,e^- 

.0  =  (  E.e^i-  E,e’’‘‘V^e^)e*' 

Where,  E^,  fi^,  E^,  are  the  constants  of  integration,  a,  p 

are  obtained  from  Equation  (3)»  and  P2,  p^  are  the  three  roots  of 
the  following  equation. 


+  P(4a.+'rs)  -K<t+sf-s(stA.){’('to.+5)  +41}  =. o 

P4  “  -  (  &!  +  s  ) 

P5  '  ai  +  s 

Shell  input;  From  the  boiindary  conditions  shown  in  Fig, 8(b),  we  get 


(52) 

(33) 

(34) 


the  transfer  function 

(ipopc-st(s) = 2  (A  e'^XCK,-  (Az-et)^ 

+  (55) 

where, 

+  0(x-  2  e\(c(,-M^-^)e  f“' 


Tube  input:  See  Fig,8(c),  The  transfer  function  is: 
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Qpcpc-ts 

Case  C-P-C-P:  The  fundasental  eqtiatlons  are, 

§+t‘=^.(e-e.) 

}  1>XL  'bJL  ^ 

11+^  =  ^,  (a -6,) 

.  '^^+11  =^(®re)+i(%-e)tA(%-e)ta(^-e) 

The  "transfer  functions  can  be  obtained  through  the  procedure 
described  before: 


Shell  Input :  See  Pig .9(b), 

Q-Cfcp-St-Cs)  =  2  (e*'»+  feXof-P*'  «*''  * 

Tube  Input:  See  Pig. 9(c),  ^  ^  J/^4 

Gcpcp-t»  (s)  =  (  "fa” 


where 


»  Pl»  Pp*  P‘5  roots  of  Equation  (32),  a. 


(38) 


(39) 


p  are  to  be 


Pl»  i'2*  *'3 

obtained  fro«  Equations  (3),  and  p^,  p^  are  given  by  Equations  (33) 
and  (34)  respectively, 

+K-p.- 


(40) 


Cases  With  More  Than  Pour  Tube-passes 

The  above  mentioned  procedures  can  be  eq\ially  applied  to  the 
heat  exchangers  which  have  more  than  four  tube-passes. 

The  procedure  is  rather  complicated  when  the  number  of  tube-pass 
is  odd • 

Por  example,  let  us  consider  the  (2n4-l)  pass  heat  exchanger: 
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/'-3.47,  C  (purecou+erjiow) 

4  88.  C-P-C 

,  -''--S.09,  P-C,C-P 

-5  m--  -  S  I  S,  P-C-P-C  .C-P-C-P 
^v'--S.19,  n  —  » 

V-S,S3.P-C-P 
^  -  6.19.  P  Cpareporalici-^low) 

\ 

-JO  -  ■  7  96.  M  Cbo+n  ilaid  mixea  cose) 


Pig. 10  Comparison  of  static  charac¬ 
teristics  of  multipass  heat  exchang¬ 
ers  for  r  =  1,  total  a^  =  total  a  = 

2. 


Pig. 11  Frequency  response  of  1-2 
pass  heat  exchangers  for  r  =  1,  = 

1. 


The  3rd  order  characteristic  equation  is; 

For  the  types  P-C-P - 

|>’  +  i-*  f(2-n+  lU  +r^}  -  S)\  aa.}-  s  (s+  a,)(^  (a,+s)+ (2-n+  i)^l} =0  , 

which  corresponds  to  Equation  (18).  The  remainder  p  values  are  given 
by  n,  n-1  I'epeated  roots  of  p^,  p^  respectively. 

For  the  types  C-P-C— ~ — 

f + t'*{(2i»+i)a+vs}  - !«{ (a,+s) -aa,}-  s  («+  (i,+s)+ftwiu) = o 

which  corresponds  to  Equation  (25).  The  remainder  p  values  are  given 
by  tt-1,  n  repeated  roots  of  p^,  p^  respectively. 


For  the  (2n)  pass  heat  exchanger,  the  characteristic  equation 
corresponding  to  Equations  (4)  and  (32)  is  as  follows  for  both 
P-C-P-C -  ^  C-P-C-P - : 

1^  +  (a,+s)'-s(s+2i,'^T(sta,)t2-n,ai}  =  o 

The  remainder  p  values  are  given  by  n-1  repeated  roots  of  p^,  p^. 

Thus,  the  dynamic  response  of  these  multi-pass  heat  exchangers 
can  be  obtained  throu^  the  use  of  these  roots.  The  order  of  the 
equations  never  exceeds  three  so  long  as  all  tube-passes  are  Identical 
and  the  number  of  shell  pass  is  one. 

If  the  number  of  tube-pass  approaches  infinity,  this  method  becomes 
very  complicated,  but  on  the  other  hand  the  temperature  gradient  in 
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Pig. 12  Frequency  response  of  1-3,  1-4  pass  heat  ex¬ 
changers  for  shell-side  Input,  r  =  1,  total  a^  = 
total  a  =  2, 


one  tubepass  becomes  almost  zero,  therefore  we  can  consider  this 

case  as  a  crossflow  type  heat  exchanger  with  both  fluids  mixed 

normal  to  the  directions  of  flow  respectively.  The  static  character- 

(1^  ■) 

istics  of  this  case  has  already  been  obtained  ,  Although  the 
dynamic  response  of  it  is  relatively  easy  to  obtain,  since  the  dead 
time  may  not  essentially  be  equal,  the  phase  characteristics  can  not 
be  obtained.  Consequently,  we  can  use  only  the  static  response 
formulas  to  obtain  the  limiting  value  of  dc  gain  in  this  case. 

Numerical  Example  of  Frequency  Response 

Static  and  dynamic  responses  are  compared  using  the  obtained 
transfer  function  formulas  neglecting  the  heat  capacities  of  solid, 

A  set  of  the  system  parameters  used  are:  r  =  1,  total  ■=  total  a  =  2, 

In  other  words,  equal  heat  exchange  conditions  for  all  heat  exchangers 
are  assitmed. 

Fig, 10  shows  that  the  dc  gain  patterns  possess  quite  systematic 
relations  each  other.  Note  that  the  response  of  a  heat  exchanger 
with  moz^  than  two  tube-passes  lies  between  C  and  P,  Also  note  that  as 
the  number  of  tube-passes  increases,  the  effects  of  counter-and 
psorallel-flow  in  individueO.  passes  are  averaged  each  other  and  the 
response  approaches  to  the  value  when  the  number  of  tube-pass  is  infinity, 
Fig,ll  shows  the  frequency  response  of  1-2  pass  heat  exchanger. 

Note  that  the  phase  curves  of  and  G^cp-ts^^'^^  ^ 
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peak  value  and  gain  ctii^es  are  oscillatory  ror  a  cerxain  xrequency 
range.  P  and  C  denote  the  responses  of  the  pure  parallel  and  counter¬ 
flow  heat  exchangers  respectively. 

Fig.l2  shows  the  responses  of  the  1-5.  1-4  pass  heat  exchangers  for 
shell  input.  Bfote  that  the  response  curves  may  he  oscillatory  about 
0)  >  1.5. 

From  Figs .11.  12.  following  resvilts  may  be  observed  and  summarised: 

(1)  The  phase  lag  increases  with  the  number  of  tube-passes. 

(2)  In  the  very  low  frequency  range,  the  gain  approaches  to  the  value 
given  by  the  dc  gain  and  phase  lag  is  small,  which  means  the  whole 
length  of  shell-pass  and  tube-pass  are  effective. 


(3)  In  the  medium  freouenc3^  range,  input  and  output  relations  have 
a  marked  effect  on  the  frequency  response  characteristics. 

(4)  In  the  high  frequency  range,  both  the  phase  ;!.nd  the  gain  curves 
become  oscillatorj"  for  oj  >  1.5,  which  might  be  due  to  the  inter¬ 
action  between  shell-pass  and  tube-pass, 

(5)  The  phase  lag  v/ould  be  larger  and  the  gain  woxild  be  more  attenuated 
when  solid  capacities  are  added  in  a  practical  heat  exchanger, 
especially  at  high  frequency  rnnge. 

Analog  Computer  Test 

Althotj^  several  analog  circixits  for  a  ptire  parallel-and  counter¬ 
flow  heat  exchanger  have  already  been  published  author 

devised  a  simpler  method  using  ordinary  analog  computer  to  investigate 
the  dynamic  response. 

Considering  case  P-C  and  dividing  the  heat  exchanger  with  no  wall 
capacity  as  shown  in  Fig, 13 (a)  to  four  equal  lumped  systems,  gives 

=  +a,(% 

e,  +a.<% 

^+(a,+i)e,=  e. 

^4(a,-n)q  =  ®.  +a4<8i: 

+(^i+ 1)  Pg  =■  67  +  9i®, 

+  =®t  +  a.  (e, -i-e,) 

r^+(24+t)fii,=®,  +a.(6!.  +  ^; 
r  ^  4(2a  + 1  t  dO 

,  T  (2a  +  f  )®^  =®,  +  a  (e^  t  e,) 

The  analog  computer  circuit  using  the  same  system  parameters  as  in 
the  former  numerical  example  is  shown  in  Fig, 13(b). 
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(a)  Frequency  responses  for  shell-  (b)  Frequency  responses  for  tube- 
side  Input.  side  Input. 

Fig. 15  Ebcperlmental  results. 


Pig. 13(c)  shows  the  frequency  response  and  the  transient  response 
for  shell-side  input.  In  some  frequency  ranges  the  phase  curve  shows 
evidently  a  peak  value,  which  would  approach  the  one  for  the  exact 
distributed  system  as  the  number  of  lumps  increases. 

Experimental  Results 

The  frequency  response  of  the  multi-pass  heat  exchanger  is  quite 
d'ifferent  than  the  others  for  1—2  pass  as  has  been  demonstrated  by  the 
numerical  example  of  Pig.ll  and  Fig,l2, 

The  author  has  made  experiments  on  a  two-pass  heat  exchanger  as 
shown  in  Fig.l4.  The  length  of  the  heat  exchanger  was  chosen  10  meters 
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closed-loop  transient  responses  to 
step-function  disturbances,  upper  diagram 
shows  the  responses  for  (a)  In  Pig. 15(a), 
lower  diagram  shows  the  responses  for  (b) 
in  Pig. 15(a). 
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Pig. 16  Indie lal  responses. 


80  as  to  cover  theoretical  peak  of  tt^  phase  shift  curve.  The  heat 
exchange  process  is  water-to~water.  The  shell  material  was  "polynite" 
in  order  to  decrease  the  effects  of  side  capacity:  tube-pass  was  thin 
(  0,3  mn  )  copper  tube.  The  test  conditions  are: 

I  for  I  for 


0,368 

8,815X10"^  kcal/m'*C 
122.7X10”^  kcal/m’C 


shell-side  input  )  tube-side  input 


10  m 


1  mn 


V  vlO  m/an,  r  *  1 
0.102  kcal/n'C 
0.682  kcal/a*C 


5.01 

0.45 

5.94 

0.92 

0.046 

70.8 

68.7 

0.25 


2.86 

0.43 

6.12 

0.81 

0.04 

62.5 

70.8 

0.22 


Tsaperature  was  recorded  by  Yaaatake-Honeywell  Electronik  Potentio- 
aeter  with  bare  Cu-constantan  wire  (  0.2  am  ^  )  as  the  temperattire 
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senslB^  elftiients,  the  respoxise  speed  of  which  is  siifficiently  high 
for  the  dyiianic  response  measurements. 

The  experimental  resxilts  which  involve  the  effects  of  series  and 
side  capacities  are  shown  in  Figs •15*  16. 

Fig.l5(a)  shows  the  frequency  responses  for  shell-side  input,  (b) 
for  tube-side  input,  Favoi-ablc  agreements  are  seen  between  the 
theoretical  and  the  experimental  restilts.  Four  cases  in  Fig, 15  each 
corresponds  to  the  four  cases  in  Fig,16  which  shows  the  transient 
responses. 

lote  that  the  phase  shift  characteristics  of  case  @  is  q\iite 
similar  to  that  of  case  ,  but  the  gain  characteristics  is  different, 
this  may  be  due  to  the  difference  between  the  heat  capacities  of  the 
input  side  fluid  and  output  side  fluid. 

Althou^  the  same  remarks  can  be  said  for  the  cases  0  >  ©  >  the 
effects  of  large  dead  time  are  remarkable. 

Control  of  1-2  Pass  Heat  Exchanger 

As  can  be  seen  from  the  experimental  resxtlts,  the  choice  of  input 
and  output  are  especially  important  for  control  applications.  This 
fact  can  be  evidently  seen  in  the  closed  loop  response  as  follows: 

Fig.lT  shows  the  responses  for  shell-side  input  when  P  and  PI 
control  actions  were  used. 

Above  figure  in  Fig .17  resxilting  from  the  case  0  in  Flg.l5  shows 
very  stable  response;  lower  one  from  the  case  0  in  Fig. 15  shows  far 
unstable  response. 
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Appendix 

A  technique  for  obtaining  the  roots  of  the  characteristic  equation 
is  to  be  described.  In  general,  the  third-order  characteristic 
equations  with  s  »=  jo)  may  be  written  in  the  form 

p^  +  p^  (  +  jn^  )  +  p  (  *2  ^“2  )  “3  ^“3  “  °  ^^5) 

where  m's  and  n^s  are  real  numbers,  and  are  ftmction  of  the  variable 

0). 

Substituting  p  «  u  +  jv  into  the  above  equation,  and  equating 

both  the  real  and  the  imaginary  parts  to  zero  yield  the  following 

two  equations. 

u^+  u^B^+  u(  Bg  -2vn^  -3v^  )  -v^b^  -▼02  +  *5  *  0 

v^+  r^nj^+  ▼(-  B2-2uBj^  -3u^)  -u^n^^  -laig  -  *  0  (46) 

Three  pairs  of  curves  can  be  drawn  on  the  u-v  complex  plane  from 
these  two  equations,  and  the  three  roots  p^^,  P2  and  p^  can  be 
found  from  the  Intersecting  points  of  the  ciuves. 

As  an  example,  consider  Equation  (4)  with  r  »  1,  and  ■  a  *  1. 

By  letting  s  the  equation  becomes 

p^  +  p^C  2  +  Jw  )  +  p(  -  1  -2^<u  )  +  Aiip-  +  ja)(w^-3  )«0  (47) 
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Flg.18  Auxiliary  curves  to  obtain  the  roots  of 
Equation  (47). 


Pig. 19  Three  roots  of  characteristic 
equation  on  complex  plane  for  r  =  1, 

=  a  =  1. 


In  order  to  find  the  values  of  p^s  for  w  *  2,  first,  the  pair  of 
equations  corresponding  to  Equation  (46)  are  found, 

u^+  2u^  +  u(  -3v^  -4v  +3  )  -2v^  +  4v  +  16  =  0  (48) 

v^+  2v^  +  v(  -3u^  -4u  -3  )  -2u^  +  4u  -  2  «  0  (49) 

By  letting  u  =  3, 2, 1,0,  and  etc,  in  Equation  (48),  two  values  of 
v*s  that  satisfy  the  equation  are  evaluated  for  each  value  of  u. 

The  result  is  the  u-curves  shown  in  solid  lines  in  Fig, 18,  The 
v-curves  are  obtained  in  likewise  manner,  and  are  shown  in  dotted 
lines  in  the  same  figure.  The  three  intersecting  points  are  the 
roots  of  the  Equation  (47), 

Althou^  the  accuracy  of  the  p’s  may  be  improved  by  expanding 
the  u,  V  scale  about  the  points  of  intersection,  the  following 
technique  used  by  the  author  is  recommended  for  its  simplicity. 

Considering  the  regeon  u  >  0,  v^'O  only,  and  expanding  scales,  one 
obtain  p^  *  0,89  +  J  2,16,  Fext,  divide  Equation  (47)  with  oj  «  2 
by  p-  0,89-  j  2,16, 
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(a) 


p^+  p(2. 89+04.16)  -3.4135+05.9448  - 

p-0.89-32.16  )  p^+  p^(2  +j2  )+  p(3  )+l6  +j2 

^  p^-  p^(0.89+j2.16) 

p^(2.89+j4.16)+  p(3  ^^^4  ) 

p^(2.89+j4.16)+  p(6.4135-09.9448) 

p(-3. 4135+ j5. 9448; 
p(-3 . 4135+ j5.9448)+15. 8787 83+02.082288 

0.121217- jO. 082288 
- (b) 

In  this  division,  0,89  and  2,16  are  modified  individually  so  as 
to  reduce  the  rresidvial  (b)  as  small  as  possible.  The  final  results  are 
0,8950  and  2,1630,  Then  from  the  quotient  (a),  obtained  through  this 
proced\ire,  the  rest  of  p  values  can  be  easily  obtained.  The  values 
s^re: 

P2  *  -  2,5320  -  i  2,0902 
P5  «  -  0,3629  -  i  2,0728 

Although  this  division  may  appear  to  be  rather  troublesome,  this 
is  more  efficient  than  obtaining  each  intersecting  points  separately. 
Applying  the  same  procedure  for  each  other  p  values  can  be 
obt6d.ned  as  shown  in  Fig,19 
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